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Abstract 

One of the primary scientific targets of current and future CMB polarization experiments is the 
search for a stochastic background of gravity waves in the early universe. As instrumental sensitivity 
improves, the limiting factor will eventually be B-mode power generated by gravitational lensing, 
which can be removed through use of so-called "delensing" algorithms. We forecast prospects for 
delensing using lensing maps which are obtained externally to CMB polarization: either from large- 
scale structure observations, or from high-resolution maps of CMB temperature. We conclude that 
the forecasts in either case are not encouraging, and that significantly delensing large-scale CMB 
polarization requires high-resolution polarization maps with sufficient sensitivity to measure the 
lensing B-mode. We also present a simple formalism for including delensing in CMB forecasts which 
is computationally fast and agrees well with Monte Carlos. 



1 Introduction 



The last two decades of cosmological observations have resulted in a model which is consistent with 
a wide variety of datasets (e.g. [H [21 El HI E] ) to an impressive level of precision. The observations 
are consistent with a flat ACDM expansion history, and Gaussian adiabatic initial conditions with 
a power spectrum which is slightly redder than scale-invariant (n s = 0.963 ± 0.012). 

A major observational frontier in the near future will be the search for a stochastic background 
of gravity waves on cosmological scales, parameterized by the tensor-to-scalar ratio r. In single- 
field slow-roll inflation, r is related to the total field variation of the inflaton in Planck units 
(Ac/) = (Mp//8-y/7r) J r l / 2 dN); "large- field" models with Planck-scale field excursions correspond 
to r > 10 -2 , whereas "small-field" models can have arbitrarily small values of r [6]. The simplest 
models of inflation, such as a power-law slow-roll potential, have r ~ 10~ 2 , and it has been argued 
that obtaining r smaller than this requires fine-tuning [7J. However, there are also models with a 
microphysical reason for requiring small r, such as the string-inspired KKLMMT model [8|, so this 
is not a firm prediction. Increasingly precise constraints on r in the next few years will sharpen this 
picture considerably. 

On the observational front, current upper limits are r < 0.36 (95% C.L.) from large-scale CMB 
temperature fl], or r < 0.72 (95% C.L.) from large-scale CMB B-mode polarization [9]. Because 
the large-scale CMB temperature contains large primary contributions (i.e. Sachs- Wolfe and Inte- 
grated Sachs- Wolfe anisotropy), the temperature constraint on r is limited by sample variance (and 
uncertainty in other cosmological parameters) and cannot be qualitatively improved. In contrast, 
large-scale B-modes are only generated by gravity waves (or other non-scalar sources, in models which 
contain them) or nonlinear effects \10\ lllj. and the polarization constraint on r is currently noise- 
limited. As instrumental sensitivity improves, the statistical error will vary roughly as er(r) oc 1/T, 
where T is the sensitivity of the experiment in detector-hours. Since detector sensitivities and array 
sizes of CMB polarization experiments are currently in a state of rapid development, polarization ex- 
periments will give the best constraints on r in the near future, and should soon probe the r ~ 10~ 3 
to 10~ 2 range where qualitative classes of early universe models can be ruled out. 

In the limit of small instrumental noise, the largest "guaranteed" source of B-mode power on large 
angular scales is gravitational lensing by large-scale structure, which converts the E-modes generated 
during recombination to a mixture of E and B |12| [T3] . The lensed B-mode power spectrum Cf lcu 
is roughly constant on large angular scales, and can thus be interpreted as an extra source of white 
noise, whose amplitude is found to be 4.4 ^K-arcmin for the fiducial cosmology used in this paper. 
When instrumental noise levels fall below this level, expriments will be lensing-limited: the limiting 
factor in constraining r will be the lensing B-mode, which acts as the dominant source of noise. 

"Delensing" algorithms have been proposed \V4\ [TC] which statistically separate the gravity 
wave and lensing B-mode signals, offering the prospect for reducing the effective noise level below 
the 4.4 /xK-arcmin floor due to gravitational lensing. These algorithms are based on the idea that, 
given a noisy observation of the CMB E-mode polarization and an estimate for the CMB lensing 
deflection field, one can form an estimate for the large-scale lensing B-mode and subtract it from 
the observed B-mode to reduce the level of lensing contamination. The required estimate for the 
deflection field is usually obtained "internally" , using a lens reconstruction estimator \17\ [TBI \19\ [20] 
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which estimates the lensing deflection field directly from observations of the lensed CMB polarization. 
High angular resolution is required for delensing (even if one is only interested in delensing B-modes 
on large angular scales), since the large-scale lensing B-mode receives contributions from E-mode and 
deflection power on small scales. In this way, after B-mode experiments become lensing-limited, the 
small and large scales will be intimately linked: the effective noise level for the large-scale B-mode 
will be determined by the fidelity with which the small-scale polarization is measured, and used to 
estimate the deflection field. 

The purpose of this paper is twofold: first, we present an approximate scheme for doing forecasts 
which include delensing. This scheme is easy to implement, computationally fast (a few CPU- 
seconds are needed per forecast) and we find that it agrees well with Monte Carlo based results in 
the literature; second, we investigate whether delensing can be made more effective by making use of 
lensing maps which are obtained externally from CMB polarization. Small-scale CMB temperature 
maps will soon be available (e.g. from the Planck satellite [21], or from ground-based experiments 
such as ACT [22J or SPT [23]) which can be used to reconstruct the CMB lensing map with reasonable 
signal-to-noise. Approximate CMB lensing maps can also be obtained using observations of large- 
scale structure (e.g. cosmic shear), but such maps will not have precisely the same redshift weighting 
as the CMB, where the contribution from large-scale structure at high redshifts is non-negligible. 
Do either of these sources of lensing information suffice (or help significantly) to delens large-scale 
CMB polarization? 

Throughout this paper, we use the WMAP7+BAO+Ho cosmology from |24j . with parameters 
n b h 2 = 0.227, Q c h 2 = 0.111, h = 0.714, n s = 0.969, A^(/s ) = 2.38 xl0~ 9 at wavenumber k = 0.002 
Mpc -1 , and r = 0.086. All power spectra were calculated using CAMB [25], using HALOFIT [26] 
to model nonlinear evolution. 

2 Preliminaries 

Gravitational lensing by large-scale structure preserves surface brightness, and therefore remaps the 
temperature anisotropy and polarization of the CMB (for recent reviews see [271I28|). The remapping 
is described mathematically by introducing a vector field d(n) (the deflection field) such that lensed 
and unlensed CMB temperature and polarization fields are related by: 



This equation is exact in the flat-sky limit, but schematic on the curved sky: here the remapping 
n + d represents an angular displacement |d| from n along a geodesic in the direction given by d, 
and there is an implicit parallel transport of the Q ±iU pseudo- vector |29j . 

To lowest order in perturbation theory, the deflection field d(n) is the gradient of a scalar lensing 
potential (i.e. d(n) = V^(n)) which can be written as a line-of-sight integral: 
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where H(z) is the Hubble factor, *&(z,x) is the Newtonian potential, D{z) denotes the comoving 
distance to redshift z, and VLk = has been assumed. In the Limber approximation, which we use 
throughout this paper, the angular power spectrum of the lensing potential is given by 

c ** - 4 f mk? p ^ k - " D(z)) (m -oh,)'- (4) 

where P<$(z, k) is the power spectrum of the potential at redshift z. The lensed B-mode is 
given in terms of the unlensed E-mode Ei 2m2 at first order in the lensing potential <j>i m by: 
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Here, the coupling coefficient ff f ^ is given by 
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A short calculation now shows that the power spectrum of the lensed B-mode is given by: 
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In writing down the form ^ of the lensed B-mode, we have made the so-called gradient approxi- 
mation, i.e. keeping only the lowest order term in the lensing potential, but this has been shown to 
be an excellent approximation for B-modes with t < 2000 |30| . 

To delens the observed (lensed) B-modes, two pieces of data are required: (1) a measurement E^ 
of the E-mode with noise power spectrum Nf E , and (2) a measurement 09£f of the lensing potential 
<p with noise power spectrum Ng. Most commonly, the measurement of the lensing potential is 
obtained "internally" from the CMB, by applying a lens reconstruction estimator, but the discussion 
in this section is general and would also apply to an "external" measurement (e.g. estimating <j) from 
observations of large-scale structure). One feature of delensing that we wish to emphasize is that 
both measurements must go from intermediate to small angular scales (roughly 30 < I < 1000), even 
if B-mode delensing is only required on large scales, since small-scale lenses contribute large-scale 
B-mode power. This is shown directly in Fig. [TJ 

Intuitively, the idea of delensing is that a lensing-limited measurement S^^m of the large-scale B- 
mode can be improved by subtracting a quadratic combination of and 4>°^ which approximates 

the lensed B-mode (Eq. (|5|) as closely as possible. More formally, the delensing estimator is given 

by 

R dcl _ nobs _ f EB ( ^ ^ 2 ^11 ^2 E °2m 2 \ ( C £ <j)°^* \ 
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Figure 1: Large scale lensing B-mode power spectrum, interpreted as a pixel noise level in ^K-arcmin, 
for varying choices of maximum E-mode multipole i^ax or deflection multipole ^ in Eq. (pj). It 
is seen that the large-scale B-mode is generated by E-mode and deflection power on a wide range of 
angular scales (roughly 30 <£ < 1000). 



If the power spectrum of the original B-mode measurement is the sum of contributions from tensor 
modes, lensing, and instrumental noise: 
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then the power spectrum of the delensed B-mode will be the sum of contributions from tensors, 
noise, and residual lensing: 
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where the residual lensed B-mode power spectrum after delensing is given by: 
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and satisfies Cf rcs < Cf lcn . 

The form ^ of the delensing estimator can be understood intuitively as Wiener filtering the 
observed E and fields, and using the filtered fields to estimate a lensed B-mode, which is then 
subtracted from the observed B-mode. A formal derivation is obtained by solving for the weights 
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on the right-hand side of (|9j) which minimize the residual power spectrum rcs . Details of this 
calculation are given in Appendix |A"} 

On large angular scales {£ < 100), we find empirically that the residual lensing power spectrum 
Cj? rcs is independent of £ to an excellent approximation^] across the wide range of instrumental 
specifications considered in this paper. This simplifies the interpretation and some implementational 
details of the forecasts which will follow: the residual lensed B-mode can be interpreted as a source 
of white noise, parameterized by a single number (the pixel noise in units /uK-arcmin) which will 
depend on the instrumental specifications. If no delensing is performed (i.e. Cf ICS = Cf lcu ), the 
pixel noise level is 4.4 /iK-arcmin in the fiducial cosmology from £jlj 



3 Forecasts 



In this paper, our basic figure of merit will be the statistical error er(r) on the tensor-to-scalar ratio 
r, in the limit of no tensor modes after delensing the large-scale B-mode, given by: 
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Let 0"o(r) denote the statistical error without performing delensing; it is given by replacing C s 
Cf leri in the above expression: 

o"o(r) = 
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We have included a parameter £ m [ n which represents the largest angular scale which can be measured 
in the presence of a sky cut and large-scale foregrounds. If this scale is large enough to include the 
signal at £ < 10 from tensor modes which enter the horizon during reionization, then the error on 
r will be very sensitive to the value of i m \ n which is assumed. This is illustrated in Fig. [2j which 
shows o"o(r) as a function of noise level for a few choices of £ m i n . For an ideal all-sky survey with 
^min = 2, the gravity wave measurement gets ~90% of its signal-to-noise from reionization scales 
{£ < 10), and «50% of its signal-to- noise from the B-mode quadrupole (£ = 2) alone. However, 
these large angular scales are very difficult to measure in suborbital experiments (although not 
impossible with long-duration ballooning) and it is not clear whether the measuring the reionization 
signal will be practical in the presence of foregrounds. Although the reionization signal at I ~ 8 
has a higher signal-to- noise ratio than the recombination signal at i ~ 60, the signal-to- foreground 
ratio is expected to be higher for the recombination signal [31], assuming that a patch of sky is 
observed which is chosen to minimize foreground contamination. For these reasons, it is currently 
unclear whether the reionization or recombination signal will offer better observational prospects for 
constraining r. 



1 To be more quantiative about this, Cf T ° B is not observationally distinguishable (in the Fisher matrix sense) from 
a best-fit constant power spectrum, even with cosmic variance limited all-sky measurements for I < 100. 
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Figure 2: Statistical error <7o(r) on the tensor-to-scalar ratio for varying noise level, assuming no 
delensing has been done. The forecast is strongly dependent on whether the reionization B-mode 
signal at I < 10 is assumed measureable in the presence of sky cuts and Galactic foregrounds. 
The value of cro( r ) levels off for noise levels < 4.4 /iK-arcmin since the observations have become 
lensing-limited. 

As remarked in the previous section, the lensed B-mode power spectrum C £ Ien and delensed power 
spectrum Cf TBS are independent of I on large scales. We will also assume that the instrumental beam 
is sufficiently small that N^ B is also approximately independent of I on large scales. In this case, if 
we define an improvement factor a = ao(r)/a(r), then a has the simple form: 



where the RHS is independent of £ (and the value of £ m i n which is assumed). 

Eq. ( p~5| ) states that a(r) factors into two pieces: a factor ao(r) which depends mainly on survey 
geometry (i.e. / s k y and £ m - m ), and a factor (1/a) which represents the improvement due to delensingj^] 
This is a very convenient simplification because it separates the issue of whether the reioinzation 
B-mode signal is measureable from the issue of how well delensing improves the r constraint. In 
the following subsections, we will show forecasts for a given various observational scenarios, and 

2 We have shown this only for a "sharp" cutoff in I parameterized by a cutoff multipole £min, but it would also 
apply in the case where the number of modes per multipole continuously drops to zero near the fundamental mode, 
e.g. because ambiguous modes [321 133] have been projected out. 
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remain agnostic on the question of whether the reionization signal can be used to obtain a clean 
measurement of r. The reader should bear in mind that the statistical error cr(r) is given by cro(r) / a, 
where the value of <7o{ r ) can be read off from Fig. [2] for a given survey geometry. 



3.1 Polarization delensing 

A B-mode measurement which extends to small angular scales (£ > 200) can be used to estimate the 
lensing potential, by means of lens reconstruction estimators which cross-correlate E and B-modes 
[T71 fTHl [TUl |2"U] . By taking this noisy estimate of 4> as the input to the delensing estimator Q, 
delensing can be performed in a way which is purely internal to a CMB polarization dataset. 

We assume that the E and B-modes have been measured with instrumental noise power spectrum 
given by 
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where Ap is the pixel noise level of the experiment (used with units /uK-radian in the above equation, 
but typically quoted in //K-arcmin). 

The noise power spectrum of the lens reconstruction obtained using the EB estimator is given 

by 
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We can plug this into Eq. (12) to obtain C« rcs , and then use Eq. (15) to obtain a forecast for a. 



However, the resulting forecast is incomplete since it does not include the improvements that could 
be obtained using iterative delensing. 

In polarization, the quadratic estimators for the lensing potential and the delensed B-mode that 
have been discussed so far can be significantly improved for low noise levels using an iterative, 
likelihood-based approach [20]. The intuition behind the iterative estimators is that lensed B-mode 



power acts as a source of noise for the lens reconstruction estimator (as can be seen from Eq. (17)), so 



that the delensed B-mode (which has less power than the lensed B-mode) can be used as the input for 
a second round of delensing with improved statistical errors, and so on iteratively. One qualitative 
difference between the two estimators is that in the limit of zero instrumental noise (N BB = 0), the 
iterative estimator can achieve perfect reconstruction of the lens potential <p and perfect delensing 
of the B-mode. As the instrumental noise improves, there is no fundamental limit to the value of r 
which can be detected [E] , unlike the case of the quadratic estimator [HI [15] , although of course 
there will be some practical limit due to foregrounds and instrumental systematics. 

Based on this intuitive picture, we propose the following procedure for forecasting B-mode mea- 
surements which use iterative delensing. After computing the residual lensing B-mode power spec- 
trum Cf rcs using Eqs. (17), ( 12 ) as described above, we recompute Cf ras iteratively until convergence, 

B B 

using the value of C f rcs from the previous iteration in place of C e lon . To state this in a completely 
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formal way, our forecasting procedure is to iterate the pair of equations 
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(18) 
(19) 



to convergence, starting by taking C £ Tl 

We have arrived at this forecasting procedure via a heuristic argument, but we can test its validity 
by comparing with the results in Table I of [16], which show values of Cf TCS obtained from Monte 
Carlo simulations of an iterative delensing estimator, for a wide range of instrumental parameters. 
We find that all entries in the table agree at the ~10% level, showing that this simple heuristic 
procedure actually provides rather accurate forecasts. Given the implementational complexity and 
computational cost of the iterative delensing estimator, this forecasting procedure is one of the main 
results of this paper. Using the optimizations from Appendix [Bj iterative delensing forecasts can be 
generated in a few CPU-seconds. 




Figure 3: Forecasted improvement a = ao(r)/a(r) in the statistical error on r due to polarization 
delensing, for varying noise level and beam. In the limit of low noise and high resolution, we find no 
limit (from delensing residuals alone) to how well r can be measured. 

In Fig.[3j we show forecasts for "internal" lens reconstruction using small-scale CMB polarization, 
for varying noise level and beam and taking ^ max = 4000 throughout. Delensing can significantly 
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improve o~(r) if the instrumental noise is < 4.4 juK-arcmin (so that the large-scale B-mode is lensing- 
limited rather than noise limited) and the beam is small (so that the small-scale lensing B-modes 
needed for lens reconstruction are measured). As the instrumental noise and resolution improve, we 
do not find any fundamental limit (due to delensing residuals alone) to how well r can be measured, 
in agreement with |16| . 

3.2 Temperature delensing 

Lens reconstruction estimators can also be used to estimate the lensing potential from CMB temper- 
ature measurements which extend to small angular scales [35 | [36 1. f39 ^ 1^0 ^ I^T ^ l^2l l^3| t^l l^5| . 
Lens reconstruction from temperature has been applied to detect CMB lensing in WMAP, in cross- 
correlation with large-scale structure (SI 07] and is expected to be a powerful source of cosmological 
information in the near future gHl HS1 [50l EH E21 [531 EH ESI [561 EZl EHl El] - It is therefore natural 
to ask whether the lens reconstruction from temperature datasets is sufficient to delens large-scale 
CMB polarization. Such a question would be relevant in a scenario where a lensing-limited B-mode 
measurement has been made on large angular scales, using an instrumental beam which is too large 
to observe the small-scale lensing B-mode needed for "internal" delensing. Note that, even if an 
"external" estimate of (ft is available from CMB temperature, polarization delensing still requires 
a measurement of E-mode polarization on intermediate to small scales (since the delensing estima- 
tor ([9]) combines E and on a wide range of scales, in order to to delens B on large scales). Therefore, 
it is not possible to avoid the use of small-scale polarization completely, but the small-scale E-mode 
measurement is much easier than the small-scale B-mode measurement needed for polarization de- 
lensing. (One possible source of the small-scale E-mode measurement would be the Planck dataset, 
which will release small-scale E-mode measurements on the full sky in the near future.) 

Since CMB lens reconstruction estimators are based on the presence of small deviations from 
Gaussian statistics induced by lensing, any non-Gaussian contribution to the CMB temperature can 
potentially bias lens reconstruction. Although a detailed characterization of such biases is largely 
unknown territory, some studies have already found significant biases from non- Gaussian signals such 
as the kinetic SZ effect |60j . Therefore, the limiting factor in temperature lens reconstruction is likely 
to be the presence of non-Gaussian secondary anisotropies (which become increasingly important as 
t increases), rather than instrumental sensitivity or resolution. However, at the time of this writing 
it is unclear what range of scales will be "sufficiently Gaussian" to use for lens reconstruction in 
practice. We will model this unclear situation in a rough way by introducing a cutoff multipole 
^max> an d assuming that temperature multipoles t < ^ ax can be used for lens reconstruction with 
the full statistical power of a Gaussian field (i.e. without introducing extra systematic error from 
secondary anisotropies), but multipoles I > ^ax are n °t used. We will furthermore assume that the 
measurements for I < ^ H have been measured with negligible instrumental noise. 

Under these assumptions, the noise power spectrum of the lens reconstruction obtained using 
the TT estimator is given by |34| : 
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Figure 4: Forecasted improvement a = cro(r)/a(r) in the statistical error on r, assuming a lensing- 
limited large-scale B-mode which is delensed using a lens reconstruction from small-scale temperature. 
We assume that the temperature measurements are cosmic-variance limited up to maximum mul- 
tipole ^nax (beyond which non- Gaussian secondary anisotropy inhibits lens reconstruction). Large- 
scale B-mode delensing also requires a measurement of the small-scale E-mode; we consider either 
the E-mode measurement expected from Planck, a future E-mode measurement with Ap = 30 
/^K-arcmin and zero beam, or a perfect E-mode measurement. 



To obtain a delensing forecast, we must also specify two pieces of data: 

1. The noise power spectrum N BB on large scales (needed to forecast a in Eq. (15)), 



2. The noise power spectrum N EE on intermediate to small scales (needed to compute C BlCB in 
Eq. pi)). 



For the first of these, we will take N BB = in order to assume that the large-scale B-mode is purely 
lensing-limited. For the second, we will consider three possibilities: either (a) Planck noise levela^J 
(b) E-mode measurements with Ap = 30 ^K-arcmin and zero beam, and (c) E-mode measurements 
with zero noise. Option (c) is artificial, since a zero-noise measurement of the small-scale E-modes 
would be accompanied by a zero-noise measurement of the lensing B-modes and one could just do 



3 When making forecasts which include Planck E-mode measurements, we use the parameters from [61) : we assume 
that the 143 GHz channel (Ap = 78 /iK-arcmin, #fwhm = 7.1 arcmin) and the 217 GHz channel (Ap = 135 /iK-arcmin, 
Afwhm = 5 arcmin) can be used with full sensitivity, while channels at lower or higher frequency are "consumed" by 
foreground cleaning. 
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polarization delensing, but we include it for the sake of having an upper bound on the possible 
improvement from temperature delensing, by neglecting uncertainty in the small-scale E-mode. 
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Figure 5: Reconstruction noise power spectra resulting from (a) temperature lens reconstruc- 
tion assuming noise-free measurements to maximum multipole ^ = 2500, (b) polarization lens 
reconstruction with Ap = 2 ^K-arcmin and #fwhm = 5 arcmin. The signal power spectrum Cf is 
shown for comparison. 



In Fig. [4j we show the improvement in a(r) that can be obtained from temperature delensing, for 
varying ^ ax and small-scale E-mode noise. It is seen that the improvement in a(r) due to delensing is 
always modest, if the small-scale CMB temperature is the source of the lens reconstruction. Obtain- 
ing an improvement larger than ~10% would require both measurements of the small-scale E-mode 
which are less noisy than Planck, and making the optimistic assumption that temperature multipoles 
out to ^ ax ~ 3000 can be cleaned of secondary anisotropy at a level where they can be used for 
lens reconstruction. This negative result can be understood in a more qualitative way by compar- 
ing the noise power spectra Ng obtained from temperature and polarization lens reconstruction 
(Fig. [5]). The lensing potential on large scales (£ < 100) is reconstructed with high signal-to- noise 
in both cases, but reconstructing lenses on small scales [l ~ 1000) requires polarization. Since 
small-scale lenses do generate B-mode power on large scales (by coupling to small-scale E- modes), 
the temperature reconstruction is of limited effectiveness in delensing the large-scale B-mode. 
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3.3 Delensing using large-scale structure 

Finally, one could ask whether an "external" estimate of obtained from large-scale structure 
observations can be used to delens large-scale B-mode polarization. One could imagine using different 
flavors of large-scale structure data (for example, cosmic shear [62J or 21-cm temperature |63 | 164 4 165]). 
with redshift and I weighting optimized to approximate the CMB lensing potential <j) as closely as 
possible. The problem encountered in the previous subsection, namely that the lens reconstruction 
from CMB temperature does not extend to small scales, is not generally an issue for large-scale 
structure, which can go to even smaller scales than the CMB polarization reconstruction. However, 
a new problem arises: the CMB lensing potential receives contributions from matter fluctuations 
at high redshift, and it is not possible to include these contributions if <j) is being estimated from 
large-scale structure data at low z. 

Let us make the crude assumption that the matter fluctuations have been observed perfectly for 
redshifts z < z max , and have not been observed at all for z > z max - (In a real survey, there will be 
a smooth transition between redshifts at which the matter density field is signal-dominated versus 
noise-dominated, but we will use a sharp cutoff to get a rough idea of what to expect.) 



A small modification to Eq. (12) for forecasting Cf Tea is needed in the large-scale structure case. 



In deriving Eq. (12), we have assumed that the measurement 0»~f of the lensing potential is the sum 



iobs 



of the true potential and an uncorrelated noise term (i.e. that (4>£m — <pim) and 4>£ m are uncorrelated 
fields). In the large-scale structure case, we write the true potential as the sum of uncorrelated 
contributions from low and high redshifts: <p£ m = 4>\° m + (fim- We then assume that <^ has been 
observed perfectly, that $f m is unobserved, and that we have an observation of the small- 

scale E-mode with noise power spectrum Np E . Given this setup, the minimum- variance delensing 
estimator and residual lensing B-mode power spectrum are given by: 

riEE pobs* \ 

R del _ R obs f EB I °l 2 ^^gmg I Ao* /o-i\ 

^imi - ^imi - 2^ Jlxlil \ n EE j_ ]\jEE I Vim 1 ~ ; ' 

/~<B TCS 




I fEB 1 2 r<EEri4>4> 
^PTZ^WW 1 U i ~ I C EE + N EE I °£ • y ZZ > 



This is shown in Appendix [Aj although the estimator (21) is easy to guess: we simply subtract an 
estimate for the lensed B-mode which is obtained by combining the Wiener-filtered E-mode and 
the observed part of the lensing potential. Note that the power spectrum Cf l ° can be calculated in 
the Limber approximation by simply restricting the line-of-sight integral Q to the redshift range 

In Fig. [6j we show the improvement in <r(r) that can be obtained from large-scale structure 
delensing, for varying z max . As in the temperature case from the previous subsection, this type 
of delensing also requires a measurement of the small-scale E-mode polarization, and we consider 
the same three cases: a Planck E-mode measurement, a measurement with Ap = 30 ^K-arcmin, 
and a zero- noise measurement. It is seen that the improvement in <r(r) due to large-scale structure 
delensing is modest, at least for survey parameters which are likely to be attainable for the next few 
generations of CMB polarization experiments. An ambitious galaxy/quasar survey which is deep 
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Figure 6: Forecasted improvement a = cro(r)/cr(r) in the statistical error on r, assuming a lensing- 
limited large-scale B-mode which is delensed using perfect large-scale structure measurements out to 
redshift z max . Large-scale B-mode delensing also requires a measurement of the small-scale E-mode; 
we consider either the E-mode measurement expected from Planck, a future E-mode measurement 
with Ap = 30 ^K-arcmin and zero beam, or a perfect E-mode measurement. 

enough to probe structure to z ~ 3 with high signal-to-noise can delens by a factor ~2, in combination 
with small-scale E-mode measurements which are a few times more sensitive than Planck. 

Qualitatively similar results were obtained in [62J, in a detailed analysis which concentrated on 
cosmic shear as the large-scale structure observable. Here, we have shown that these conclusions are 
generic to any observable which does not probe the matter distributions at high redshift. Even with 
perfect knowlege of the matter fluctuations out to redshift ~3, only relatively modest improvements 
in a(r) can be obtained. 

4 Discussion 

In this paper, we have presented forecasts for delensing the gravity wave B-mode, aimed at assessing 
whether high-sensitivity measurements of CMB polarization on small scales are a requirement, or 
whether the lensing map can come from a different source. We have also presented a simple frame- 
work for including delensing in forecasts, which is approximate but agrees very well with Monte 
Carlo simulations. 

One technical but important point is that calculations throughout this paper have treated all 
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fields as Gaussian, even "composite" fields which are formed from products of others (such as the 
lensed B-mode, which is a composite of E and <p). In more detail, we have made three approximations 
as follows. First, the Fisher forecasts for a(r) and <Jo(r) (Eqs. ( [131 ), (14)) have assumed that the 
bandpower covariance of the lensed B-mode B^ and delensed B-mode BIT! are related to the power 
spectrum in the same way as for a Gaussian field. For the lensed B-mode, it has been shown [66J that 
this is a good approximation for narrow bandpowers, but breaks down for bandpowers wider than the 
coherence scale of the lenses (A£ > 200). Since the effective bandpower for constraining r is narrow 
(A^ ~ 40), the approximation should be accurate in this case. It seems plausible that the same 
qualitative statements are true for the delensed B-mode, but we defer investigation of this issue to 
further work. The second approximation appears in Eq. (12): when calculating the residual B-mode 
power spectrum C^ Tes , we have used Wick's theorem to compute expectation values of products of the 
fields B un \ E ohs , cp ohs . More precisely, we have neglected the (B unl E ohs (j) ohs ) three-point function and 
the ^E ohs E ohs (j) ohs (j) ohs ) connected four-point function. The validity of this approximation is unclear 
in the case where cp ohs is obtained by applying a quadratic lens reconstruction estimator to the CMB; 
in this case the 3-point and 4-point terms which have been neglected will "expand" to a nontrivial 
CMB 4-point or 6-point function. The situation becomes even murkier if the delensing is iterative, 
since iteration will lead to Appoint functions of even higher order. The third approximation we have 
made is treating the lensing potential <p as a Gaussian field. We have included nonlinear corrections 
to the power spectrum Cf , but nonlinear evolution also generates higher A-point functions which 
have been neglected. A complete study of these approximations and potential biases is outside the 
scope of this paper, and will be pursued separately [67]. As remarked in £3.1 the close agreement 
between the forecasts in this paper and the Monte Carlo results from [16J suggests that the impact 
of these approximations is small in forecasts, but may be important in actual data analysis. 

We considered the possibility of using the CMB lens reconstruction from small scale temperature 
to help delens large-scale polarization. We found that this is not a promising strategy and that it 
will be difficult to achieve a reduction in a(r) which is better than ~ 10% with realistic assumptions 
about small-scale secondary temperature anisotropy. The limiting factor is poor signal-to-noise in 
the TT lens reconstruction on small angular scales (£ ~ 1000). 

We also considered the possibility of using large-scale structure to help delens the CMB. While 
the forecasts here are not especially encouraging, they are not as grim as the temperature delensing 
case. At some point in the future, it may be interesting to obtain a modest (say 50%) improvement 
in cr(r) using large-scale structure, as a proof of concept of the delensing idea. (Large-scale structure 
delensing may also be easier in practice than polarization delensing, since there are fewer systematics 
which will cross-correlate the two datasets.) Going beyond this would require small-scale polarization 
measurements and polarization delensing. Another futuristic possibility is a 21-cm experiment with 
sufficient sensitivity and resolution to reconstruct the lenses out to redshift ~10 |63l l64"| 168] , but 
this is unlikely to be available for the next few generations of CMB polarization experiments. We 
conclude that polarization delensing is the only prospect in the forseeable future for probing the 
gravity wave signal significantly beyond the noise floor imposed by gravitational lensing. 
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A Appendix: minimum-variance estimators 



In this appendix, we derive the minimum-variance delensing estimator and the residual B-mode 
power spectrum, in the lens reconstruction case (Eqs. (12)) and the large-scale structure case 
(Eqs. @, d22b). 
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Considering first the lens reconstruction case, we assume that the observed B-mode is the sum 
of contributions from tensor modes, noise, and lensing: 

nobs otens . jsnoisc , Rlcn 

a t 1 m 1 — D lvm\ n hmi + D t 1 m 1 

_ R tens | rjnoise , \ ^ fEB ( ^2 ^ \ rp* A* /'OQA 

' \ mi m2 m / 

l 2 m 2 tm \ / 

and that we have noisy observations Ef^^ 0°^ s of the small-scale E-mode and lensing potential, 
with noise power spectra N EE and Nf^. We consider a general delensing estimator of the form 



i 2 m 2 ira \ 
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with weights ge^e to be determined. The form of the second term (in which g does not depend 
on m, mi, 7712) is the most general form which is bilinear in E,cp and satisfies rotation invariance. 
A short calculation shows that the power spectrum of the delensed B-mode is given by C £ dcl = 

C B tcns + N BB + C B ves ^ where 

riBrcs _ 1 \ " \( f EB f EB* f EB * „ f EB*\ r ,EE n <j>(j> 

1 e a t 
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Now we determine gi x t 2 i by minimizing the residual power spectrum Cf res . Differentiating (25) with 
respect to g, we find that the minimum occurs at 

/ c ee \ / C H \ 
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and the value of Cf ICS at minimum is given by 
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(27) 



This completes the derivation of Eqs. Q, (12) and concludes the lens reconstruction case. 

Moving on to the large-scale structure case, we assume that the lensing potential is a sum of 
uncorrelated terms 

<t>im = 4>tn + <& , (28) 

where (fi^ m has been observed without noise, and that we have a noisy observation E'^ n of the 
small-scale E-mode with noise power spectrum N EE . 
We consider the delensing estimator: 
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In this case, the residual B-mode power spectrum is given by 

,B 
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Solving for the weights hg^i which minimize C EieB , we get 
and residual B-mode power spectrum 



(30) 
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(32) 



This completes the derivation of Eqs. (21), (22) and finishes the large-scale structure case. 



B Appendix: fast real-space expressions 

The forecasts in this paper require repeated calculations of the lensed B-mode power spectrum 



(Eq. Q) and lens reconstruction power spectra in the TT (Eq. 20) and EB (Eq. (17)) cases. We 



have written expressions for these quantities as harmonic-space sums whose computational cost is 
In this appendix, we will present equivalent real-space expressions which can be evaluated 
exactly with cost O(^max)- The method to convert these expressions to real space was originally 
proposed in |70j . This optimization is very convenient since forecasts can be done in a few CPU- 
seconds for a given set of instrumental specifications, permitting rapid exploration of parameter 
space. 

We first note that the F symbol defined previously in Eq. ([7]) can be equivalently written 



r £xl2£z 



(2*i + 1)(2* 2 + 1)(4)(4 + 1)(2£ 3 + 1) 



167T 

V(.h-s)(£ 2 + s + l) ( 
+V(t2 + s)(£ 2 -s + l) 



S 8 + 1 -1 I 

( h h h 
[-8 8-1 1 



(33) 



We also have the identity 
■ i 
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where d \ , are Wigner d- functions such that s± + s 2 + S3 = s'i + s' 2 + s' 3 



0. 
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Armed with identities (33), (34), a long but straightforward calculation shows that Eq. Q for 
the lensed B-mode power spectrum is equivalent to: 



C 



d(cos6) (C£_3(*)d(0) + 2^(6)^(9) + Ct-i(e)d{0)) 4,-2(0) , (35) 
where the correlation functions which appear are defined by 
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(36) 
(37) 
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(39) 



To obtain a fast algorithm for computing C» lcn , we first note that the integral on the RHS of (35) 
can be evaluated exactly using Gauss-Legendre quadrature with |~(3£ max + l)/2] points, since the 
integrand is a polynomial of degree < 3£ max . We precompute the correlation functions in Eqs. (36) 
(39) at the Gauss-Legendre quadrature points with computational cost O(^max); evaluating the 
Wigner (i-functions by upward recursion in I. The integral can then be done for all values of I with 
the same C(^max) cost, using the same recursion. This completes the fast algorithm for C^ leiL . 

Analogously, a fast algorithm for the EB lens reconstruction noise can be obtained using the 
expression (equivalent to Eq. (17)): 
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For the TT lens reconstruction noise , we use the expression (equivalent to Eq. (20)): 

tt£(£ + 1) J d(co S d) (CS.WCi T iW-C?iWCo T iW)^-i,-iW 
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with correlation functions defined by 
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(45) 



(46) 
(47) 
(48) 



Fast algorithms can also be given for lens reconstruction noise power spectra which arise from TE, 
TB, and EE estimators, but we omit these cases here since they are not needed in this paper. 
Indeed, fast algorithms exist for calculating the noise power spectra for many quadratic estimators 
of statistical anisotropy, provided that the weights which they use are separable. 
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